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摘 要 : 青藏 高 原 是 全 球 气候 变化 最 敏感 的 地 区 之 一 。 计 算 青 藏 高 原生 态 系统 净 初 级 生产 力 (Net 
Primary Productivity , NPP) 对 精确 估算 全 球 碳 循环 具有 着 重要 意义 。 基 于 CEVSA 模型 ,利用 M-K 
趋势 检验 法 Sen’s 斜率 估计 法 及 Pearson 相关 系数 法 ,分 析 了 2000—2014 年 青藏 高 原生 态 系统 的 
净 初 级 生产 力 时 空 变化 特征 。 结 果 表 明 :(1) 青藏 高 原 高 寒 生 态 系统 净 初 级 生产 力 在 空间 分 布 上 


表现 出 由 东南 向 西北 减 小 的 趋势 ,在 东部 及 东南 部 的 森林 区 NPP 在 600 ~1 200 gC - m^ «a^! z 
间 , 中 部 草原 和 草 旬 区 NPP 在 200 ~400 eC mm- 
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“a 之 间 , 西 部 和 北部 荒漠 区 , 受 水 热 条 件 的 限 


fi] NPP 很 小 ,该 趋势 与 水 热 分 布 趋势 基本 一 致 。(2) NPP 年 际 变化 与 多 年 平均 气温 呈正 相关 ,与 
降水 量 呈 负 相 关 。NPP 与 气温 呈正 相关 的 地 区 面积 占 研究 区 总 面积 的 82.249% ,与 降水 量 呈 负 相 
关 的 地 区 面积 占 49.31% ,表明 气温 是 影响 植被 NPP 空间 分 布 的 主要 因子 。(3) 近 15 ak, FH 
高 原 NPP 整体 呈 增 加 趋势 ,与 气温 趋势 变化 一 致 ,降水 量 表现 出 微弱 的 减少 趋势 ,气温 的 增加 伴随 


降水 量 的 减少 是 青藏 高 原 WPP 缓慢 增加 的 主要 原 


因 。 因 此 ,准确 描述 NPP 对 气候 变化 响应 的 能 


力 将 使 我 们 能 够 深入 理解 陆地 生态 系统 应 对 全 球 变化 做 出 的 反应 。 
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自 20 世纪 30 年 代 以 来 , 随 着 经 济 的 发 展 , 人 类 
加 大 了 对 自然 的 开发 力度 ,煤炭 、 石 油 等 能 源 的 广泛 


CEVSA 模型 ; 净 初 级 生产 力 ; 气候 变化 ; 青藏 高 原 


应 是 我 们 了 解 陆地 生态 系统 碳 吸 收 能 力 和 相关 生物 
地 球 化 学 过 程 的 基础 5 。ZHAO 等 中 的 研究 揭示 


使 用 ,导致 许多 地 区 的 生态 系统 都 受到 不 同 程度 的 
干扰 ,青藏 高 原 是 全 球 最 敏感 的 区 域 之 一 。 高 山 
生态 系统 在 全 球 碳 循 环 中 具有 重要 作用 ”| ,在 一 定 
程度 上 具有 减缓 大 气 CO, 浓度 增加 和 气候 变化 的 
VERA?) 。 许 多 研究 表明 ,青藏 高 原 对 北半球 的 当地 
天 气 和 气候 以 及 大 气 环流 都 产生 了 深远 的 热力 和 动 
力学 影响 ” 。 气 候 变 化 是 影响 碳 循环 的 因素 之 一 ， 
鉴于 高 山 生 态 系统 生产 力 对 生长 季节 变 暖 的 响应 ， 
监测 和 量化 青藏 高 原 高 山 生 态 系统 中 积累 的 碳 量具 
有 十 分 重要 的 意义 "|。 

净 初 级 生产 力 (Net Primary Productivity , NPP ) 
是 指 速率 过 程 , 即 每 天 ,每 周 或 每 年 产生 的 植物 物质 
的 量 ( 净 初 级 生产 量 ) ,是 陆地 生态 系统 碳 通 量 的 关 
键 组 成 部 分 。 佑 算 NPP 并 评估 其 对 气候 变化 的 响 
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了 ,年 际 NPP 的 变化 与 大 气 中 C0, 浓度 有 明显 的 负 
相关 关系 ,表明 NPP 在 调节 大 气 CO, 浓度 中 起 着 举 
足 轻 重 的 作用 。 

因此 ,NPP 对 于 生态 系统 碳 通 量 的 汇集 、 调 节 
全 球 气候 变化 扮演 着 重要 角色 "] 。 青 藏 高 原 (QTP) 
是 地 球 上 海拔 最 高 且 最 大 的 高 原 , 被 称 为 世界 “第 
三 极 ” ,是 东亚 地 区 重要 的 水 库 ,在 调节 气候 变化 和 
水 资源 分 配方 面具 有 重要 作用 "| 。 近 年 来 , 随 着 数 
据 资 源 易 获 取 和 多 种 研究 方法 的 改进 ,国内 外 学 者 
应 用 各 种 算法 和 模型 对 青藏 高 原 植被 NPP 进行 了 
估算 。PIAO 等 中 等 利用 遥感 数据 ,通过 CASA 模型 
模拟 得 到 了 青藏 高 原 NPP 的 时 空 变化 特征 。LI 
等 (中 基于 CASA 模型 模 了 青藏 高 原 NPP, 并 以 净 
初级 生产 力 为 基础 定量 评估 了 和 气候 变化 和 人 类 活动 
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在 青 藏 高 原 荒 漠 化 过 程 中 的 相对 作用 。 以 上 研究 时 
段 主 要 集中 在 1980 年 以 后 的 20 a, 近 十 几 年 来 ,对 
NPP 变化 趋势 及 特征 分 析 往 往 集 中 在 全 国 尺 度 上 ， 
对 青藏 高 原 研究 较 少 。 因 此 ,对 这 方面 的 研究 有 必 
要 进行 补充 和 完善 。 

CEVSA 过 程 模 型 基于 生理 和 生态 过 程 构建 , 结 
合 光 合作 用 、 燕 散 发 .呼吸 作用 和 干 物质 分 配 来 估算 
NPP 能 够 在 理论 上 详细 解释 净 初 级 生产 力 "| 。 基 
于 此 ,本 文 应 用 CEVSA 模型 M-K 趋势 检验 法 、 
Sen 's 斜 率 估计 法 及 计算 Pearson 相关 系数 法 分 析 其 
时 空 分 布 格局 和 多 年 变化 趋势 。 模 拟 青藏 高 原 植被 
NPP 时 空 变化 及 其 对 气温 和 降水 量 的 响应 ,能 够 准 
确 监测 和 量化 青藏 高 原 高 山 生态 系统 积累 的 碳 量 ， 
售 算出 青藏 高 原 植被 NPP 在 气候 变化 的 影响 下 对 
全 球 碳 循 环 的 贡献 。 


1 研究 区 概况 与 方法 


1.1 研究 区 概况 

青藏 高 原 的 面积 为 2.6 x 10° km ,由 于 地 形 多 
样 ,气候 非常 复杂 。 东 南部 植被 类 型 以 森林 为 主 , 气 
候 温暖 湿润 ;西北 部 多 为 戈壁 气候 寒冷 干燥 。 在 整 
个 地 区 ,年 平均 气温 范围 为 -5.6 ~17.6 C ,日 温差 
为 14 ~17 % ,降水 量 分 布 不 均匀 。 东 南部 地 区 和 西 
北部 地 区 的 年 降水 量 分 别 超过 4 000 mm 和 17.6 
mm。 在 西风 环流 和 高 原 地 形 的 影响 下 ,年 平均 风速 
>3.0 ms ,大 风 ( 即 风速 高 于 17 m + s) 每 年 
RÆ > 50", 
1.2 CEVSA 模型 及 原理 

基于 过 程 的 生态 系统 模型 CEVSA ( Carbon Ex- 
change between Vegetation, Soil and Atomasphere ) 模 
型 ,用 于 模拟 植被 ER KEZA ADK CRUS T 
循环 及 对 环境 变化 的 响应 ”” 。 它 由 3 个 子 模型 
构成 ,分 别 为 :生物 物理 子 模型 用 于 估计 植物 一 土 
壤 一 大 气 之 间 水 热 交 换 土壤 含水 量 和 气孔 传导 等 
过 程 ;植物 生理 生长 子 模型 计算 植物 光合 作用 呼吸 
作用 、 氮 吸收 速率 以 及 碳 氮 在 植物 各 器 官 之 间 分 配 、 
积累 、 周 转 和 凋落 的 量 ; 土 壤 碳 氮 转 化 子 模 型 用 于 估 
计 土 壤 有 机 质 分 解 与 转化 和 有 机 损 矿 化 等 过 程 。 
CEVSA 模型 中 使 用 的 生物 学 和 生态 学 原理 ,方程 和 
参数 来 自 许 多 实验 室 和 现场 实验 和 观察 。 本 文 介 绍 
了 CEVSA 模型 的 主要 过 程 及 其 原理 ,关于 模型 的 详 


细 结 构 示 意图 如 下 : 
1.2.1 光合 作用 和 和气 筷 导 度 FARQUHR 4x 4 
绍 了 一 个 生物 化 学 过 程控 制 光合 作用 的 模型 。 该 模 
型 考虑 了 CO, 浓度、 辐 照 度 及 光 呼 吸 速率 对 净 光 合 
的 有 影响。 光合 速率 A, 取决 于 三 个 过 程 最 小 速率 : 
(1) REEK, W, : Wi F Rubisco 酶 的 量 、 动 
力学 性 质 及 活性 等 。 
(2) 电子 传输 速率 限制 ,WW 受 限 于 卡尔 文 循环 
二 磷酸 核 酮 糖 的 再 生 速 率 , 而 该 过 程 受 限于 电子 传 
递 速率 。 


辐射 平衡 热 传 递 
ROR 径流 


oka. m. aaun] 


图 1 CEVSA 模型 结构 示意 图 
Fig.1 Schematic diagram of the CEVSA model structure 


(3) 磷酸 丙 糖 (triose phosphate ) 利 用 率 对 光合 
速率 的 影响 ,该 限制 主要 发 生 在 高 C0, 浓度 或 高 辐 
照度 时 。 相 应 羧 化 速率 用 W, 表示 。 

每 层 叶 片 的 光合 作用 通过 考虑 土壤 水 分 对 气孔 
导 度 的 调节 作用 将 Farquhar 模型 及 改进 的 Ball-Ber- 
ry 模型 相 结 合 模 拟 得 到 ,NPP enr AR pd EE 
分 配 。 利 用 辐射 分 配方 法 ,将 叶片 的 光合 作用 上 升 
到 冠 层 尺度 ,对 叶 面 积 进行 分 层 , 层 数 等 于 冠 层 叶 面 
积 指 数 。 基 于 Beer’s Law 可 以 计算 得 到 每 一 层 叶 面 
积 指 数 La 的 平均 辐 照 度 (Trradiance, 1) ; 

T=he- (1) 
式 中 :7m 为 冠 层 和 人 射 辐 照 度 ;K 为 消光 系数 (设置 为 
0.5) ; 辐 照 度 7 用 于 计算 电子 传输 速率 。 分 配 到 冠 
层 每 一 层 叶 片 的 氮 与 叶片 平均 辐 照 度 成 比例 : 


I 
N-N, A (2) 
式 中 :NN 为 叶片 氮 吸 收 速率 ;Ni 为 总 氨 吸 收 速率 。 


1.2.2 和 氢 吸 收 和 自 养 呼吸 ” 碳 氮 之 间 的 耦合 作用 
包括 氮 对 呼吸 作用 的 限制 及 土壤 碳 氮 分 解 等 过 程 。 
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许 ” 洁 等 :2000 一 2014 年 青藏 高 原 植被 次 初级 生产 力 时 空 变 化 及 对 气候 变化 的 响应 


植物 的 氮 含 量 与 植物 的 最 大 光合 速率 有 明显 的 正 相 
关 性 , 氮 吸 收 对 气温 的 依赖 性 (NN;) 可 表达 为 : 


190N 
max ~ 306 +N (3) 
NOR TF) 
Np = Ty, -205.9 ky ( T) (4) 


1 + e(0-008317;, ) 


式 中 :4 为 最 大 光 饱 和 的 光合 速率 ;w 为 某 一 层 叶 
片 的 氮 含 量 ;Ny 则 为 植物 的 总 氮 量 吸收 量 ; ui 、w、 
us 和 kr(7T) 均 为 气温 的 响应 函数 ;7 为 绝对 气温 。 
叶片 的 氮 含 量 (N) 还 决定 了 叶片 的 维持 呼吸 
(Ra): 


N 
Rm > 50 


eO “ea; (5) 
SUP sr, (T) Ar, CT) 400 ^C AY Ws] P BRL o 

1.2.3 土壤 碳 和 氮 动 态 CEVSA 模型 参考 CENTU- 
RY 模型 来 模拟 土壤 碳 氮 分解 过 程 。CENTURY Ti 
型 模拟 了 土壤 系统 C.N.P.S 的 有 机 和 无 机 动态 变 
化 ,CEVSA 模型 仅 采 用 了 土壤 碳 部 分 。 土 壤 有 机 质 
可 分 为 地 表 凋 落 物 碳 库 地 下 凋落 物 碳 库 .微生物 碳 
P .慢性 碳 库 及 惰性 碳 库 ,通过 一 级 反应 速率 模拟 各 
碳 库 的 碳 分 解 忆 ] 。CEVSA 模型 通过 以 下 几 个 函数 
模拟 氮 对 土壤 有 效 碳 分 解 的 限制 : 

微生物 可 利用 的 潜在 碳 (Ca) 及 潜在 氮 (Na): 


Ca= 5,0MK, (6) 


Na= >,OM.K,(CN,) ^! (7) 


式 中 :i=1 ~8 5 OM 为 碳 库 ;K 为 潜在 周转 速率 (考虑 
了 气温 .土壤 水 分 等 限制 因素 ) ;CN 为 CAN 比 。 
土壤 异 养 呼 吸 (HR) 可 用 下 式 表达 : 
HR = 5,0M,K,, (1-8) (8) 


式 中 :Ki 是 实际 各 个 土壤 有 机 碳 库 的 衰变 率 ;e 为 
同化 速率 。 

1.2.4 LEKAR ARK LHESOKHORIEWOB 
RRI TOR PR i AIL SE D BR CF BY), HE AY 
ORF RAK HA RAR BR : 

SMOS,, = min | [SMOS,,_, + WATI -AET] ,SMOS,,,| (9) 


式 中 :m 表示 时 间 步 长 (月 ); WATI (mm. + month ^) 
为 总 的 水 分 输入 ;包括 降水 和 降雪 融化 ; SMOS,,, 为 
饱和 土壤 水 分 。 当 土壤 水 分 高 于 SMOS,, 时 ,水 分 作 
为 径流 损失 。 


降雪 与 降水 之 间 的 分 配 被 描述 为 : 
0.0 TEM z2 C 


SNFC - ae 2'C»TEM» -5'C (10) 


1.0 TEM x -5 C 


式 中 :SNFC 是 融雪 作为 降水 的 部 分 ;TEM 是 每 个 月 
平均 气温 。 

基于 ABER 等 .理论 ,当月 均 温 高 于 1 CHE, 
积 雪 开始 融化 (SVME ) ,最 大 融雪 量 为 当前 积 雪 量 
最 大 值 : 


SNME =45TEM (11) 
因此 ,土壤 水 分 输入 为 : 
WATI = PRCP(1-SNFC) + SNME (12) 


Federer 研究 ,水 分 损失 (4E7T) 为 供应 函数 S 及 
需求 函数 D Bp MELOS: 


SMOS 
5 2 ET, S MOS (13) 
AET = MIN(S,D) (14) 
SUP LET, ON OK ZR BOR « 
ZR S RRL D , 4 BE REDE : 
sS R, 
D= i +y | L (15) 


式 中 :RR, 为 净 辐 射 ;s 为 饱和 水 气压 随 气 温 变 化 斜率 
(Pa + K^) vy 为 干 湿 指 数 ;L 为 水 分 蒸发 潜 热 。 
1.3 输入 数据 来 源 、 参 数 获取 和 模型 运行 

CEVSA 模型 输入 数据 中 10 d 时 间 步 长 的 平均 
气温 和 降水 量 ( 旬 内 数据 累加 ) 数 据 来自 中 国 气象 
数据 网 CNOS. 1 数据 , 云 量 数据 和 相对 湿度 数据 来 
自 于 欧洲 中 期 天 气 预报 中 心 的 ERA-Interim 再 分 析 
资料 ,后 运用 MATLAB 软件 重 采 样 到 0. 1° x0. 1? 分 
BE, KA CO, 浓度 数据 来 自 地 球 二 氧化 碳 网 站 
( https ://www. co2. earth/)。 土 壤 类 型 和 质地 数据 
取 自 世界 土壤 数据 库 HWSD 中 国 数据 库 , 对 其 进行 
数字 化 并 重 采 样 到 0.1°。 植 被 类 型 分 布 数据 来 自 1 
km 分 辩 率 的 全 球 土地 覆盖 数据 集 ( European Com- 
mission, Joint Research Centre 2003. http://biova. 
Jre. Ec. europa. eu/products/g1c2000/21c2000. php) 。 

CEVSA 模型 包含 43 个 参数 和 11 个 初始 的 状 
态 变量 “' 。 在 CEVSA 模型 应 用 于 区 域 和 全 球 的 模 
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拟 时 ,初始 的 状态 变量 ,如 植被 碳 .土壤 碳 EAL 
比 、 初 始 的 土壤 含水 量 等 来 源 于 静态 模型 的 输出 。 
与 土壤 质地 有 关 的 参数 包含 田间 持 水 量 、 饱 和 含水 
Ht ES AAR .土壤 分 解 的 最 佳 持 水 量 等 参数 ,这些 
土壤 水 分 参数 由 土壤 颗粒 组 成 确定 ,依据 RAICH 
45 MCGUIRE 45" CAO 4&' 文献 中 提供 的 数 
据 进 行 估 算 。CEVSA 包含 13 种 土地 覆盖 类 型 ,与 
植被 类 型 相关 的 参数 主要 是 根 、. 蔡 、. 叶 的 周转 时 间 , 
根 ` 蔡 . 叶 的 分 配 比 例 和 比 叶 面 积 等 ,主要 来 源 于 已 
发 表 文献 的 估计 。 本 研究 选取 青藏 高 原 的 土地 和 窗 盖 
类 型 ,如 与 最 初 模型 中 已 有 土地 履 盖 类 型 直接 对 应 
的 ,就 使 用 原 有 参数 ,不 直接 对 应 的 植被 类 型 ,本 研 
究 通过 查阅 文献 和 模型 调试 的 方法 获取 相对 应 植被 
类 型 的 相关 参数 。 详 细 的 参数 列表 和 取 值 方法 可 以 
参考 已 发 表 的 文章 ,如 CAO a6 09 40 ZHANG 2%) 。 
本 研究 中 ,模型 的 空间 分 辩 率 为 0. 1?。 首 先 ， 
我 们 使 用 54 a 平均 气候 (1961 一 2000 年 ) 和 固定 的 
1960 EKA CO, 浓度 值 来 运行 模型 ,直到 碳 库 达到 
平衡 , 即 NPP 凋落 物产 生 量 (7L7) 与 异 氧 呼吸 (RH) 
相等 ,然后 使 用 1961 一 2014 年 期 间 的 气象 数据 和 大 
气 二 氧化 碳 浓度 数据 进行 模拟 。 为 了 消除 使 用 不 同 
初始 值 对 模拟 的 NPP 不 确定 性 的 影响 ,模型 利用 
1961 一 2000 年 的 数据 模拟 来 确定 2000 一 2014 年 期 
间 的 实际 初始 值 。 
1.4 时 空 变化 及 对 影响 因素 的 分 析 方 法 
1.4.1 检验 序列 变化 趋势 ”本 文 应 用 Manner-Ken- 
dall( M-K) 非 参数 检验 法 来 计算 NPP 降水 量 和 气 
温 等 随时 间 变 化 的 趋势 。 假 设 有 nn 个 样本 (x ,…， 
x.) ITEFA, XI FAH k jsn, H kAj, x, I x, 
HAER, PROSE Ss, ARF : 


p Sena; m) (16) 
+1 (x =m,) >0 

HP Sgn(x,-x,) = 0 (a,-«) =0 (17) 
-1 («,-%,) <0 


S 为 正 态 分 布 ,均值 为 0 ,方差 Var(s) =n(n - 
1) (2n 45)/18, ?4 n >10 时 ,标准 的 正 态 统计 变量 
通过 下 式 计算 : 


(18) 


对 于 统计 值 Z 来 说 , >0 时 ,表示 序列 增加 ; 反 
之 则 减少 7, 

Sen's 斜率 佑 计 法 是 选择 通过 成 对 点 的 所 有 线 
的 斜率 的 中 值 来 稳健 地 将 线 拟 合 到 平面 中 的 采样 点 
(简单 线性 回归 ) 的 方法 ,该 方法 能 很 好 的 反应 序列 
趋势 的 变化 程度 ,Sen's 斜率 天 的 公式 如 下 |: 


K =median(Q;) (19) 
nur Eu 


式 中 :i 为 常数 取 值 为 (1,2,…,N) ;N 为 时 间 序 列 长 
度 ;x; Mx, 则 分 别 是 7 入 的 序列 ; 且 满 足 7 > ko 
1.4.2 计算 相关 系数 ”由 于 本 文 NPP 是 通过 驱动 
气象 因子 在 CEVSA 模型 中 模拟 得 到 的 ,因此 ,NPP 
与 气温 和 降水 量 之 间 有 一 定 的 相关 性 。 为 深入 分 析 
NPP 对 降水 量 和 气温 的 响应 ,本 文采 用 皮尔 逊 相关 
系数 的 方法 ,计算 气温 和 降水 量 与 NPP 的 相关 性 。 


2 研究 结果 与 分 析 


2.1 基于 观测 数据 的 CEVSA 模型 估算 结果 验证 

CEVSA 模型 已 经 从 站 点 到 区 域 范 围 进行 了 验 
证 。 在 站 点 尺度 , 顾 峰 雪 等 ”使 用 不 同 植被 功能 > 
型 的 涡 动 协 方差 通 量 观测 数据 来 验证 模型 性 能 。 本 
文中 由 于 缺乏 研究 区 内 植被 NPP 的 观测 值 ,加 之 研 
究 区 内 观测 站 点 数量 少 ,对 于 模型 估算 进行 验证 是 
很 困难 的 。 本 文 仅 从 FLUXDATA 上 获取 了 研究 区 
Vj 4 ME (2004—2005 年 ) 海北 草 甸 (2002 一 2004 
年 ) 海北 湿地 (2003 一 2005 年 )3 个 观测 站 点 进行 
模拟 并 验证 。 表 1 是 3 个 站 点 观测 与 模拟 的 年 总 净 
初级 生产 力 (NPP) 的 对 比 。 

2004—2005 年 当 雄 模拟 值 与 观测 值 的 相关 性 
分 别 为 0.77 和 0.82; 2002—2004 年 海北 草 旬 模 拟 
值 与 观测 值 的 相关 性 分 别 为 0. 87、0.94 和 0. 84; 
2003 一 2005 年 海北 湿地 模拟 值 与 观测 值 的 相关 性 
分 别 为 0.83 .0.85 和 0.87。 图 2 分 别 为 当 雄 海北 
草包 海北 湿地 总 体 观 测 值 与 模拟 值 的 线性 拟 合 散 点 

。 总 体 来 看 , 3 个 站 点 中 模型 估算 的 植被 VPP 与 

观测 NPP 值 有 较 好 的 一 致 性 , 且 均 在 合理 范围 之 内 。 
2.2 区域 多 年 平均 气温 、 年 降水 量 和 年 NPP 的 时 
空 分 布 
2.2.1 Ai a 青藏 高 原 年 均 气 温 、 年 均 降 水 量 及 
植被 年 均 NPP 的 时 间 变 化 ” 近 15 a 来 青藏 高 原 
NPP 总 量 的 平均 为 0.56 Gta! (C), 是 全 国 NPP 
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R1 3 个 站 点 年 总 NPP 观测 值 与 模拟 值 的 对 比 / gC- m? «a^ 


Tab.1 Comparison of total NPP observations and simulated values of three stations / gC - m ^ - a^! 
ee 观测 值 模拟 值 
观测 站 = = 
2002 年 2003 年 2004 年 2005 年 2002 年 2003 年 2004 年 2005 年 
EE = = 405 310 i i 340 329 
海北 草 负 414 475 568 一 362 400 468 一 
海北 湿地 一 523 580 573 = 409 430 410 
7 40, (a) Sit 了 人 b) 海 北 草 旬 1.40, (c) 海 北 湿地 
T © 50 T 
= = =| 
30 上 T 40 ky 30 | 
Q, 20 3 220r 
i im 20 Es 
$10; EN RZ 10 je 
im R°=0.71 æ 10 R2-0.70 E R°=0.81 
3 3 E ale 
z 0 Lm we 。 z 0 e. zZ a 
0 10 20 30 40 50 60 0 20 40 60 80 100 120 0 20 40 60 
NPP3Wil& / gC * m? (108) NPP3Wil& / gC * m? (10d) NPP3WiÉ / gC * m? (10d) 


图 2 2004—2005 年 当 雄 2002 一 2004 年 海北 草 旬 、.2003 一 2005 年 海北 湿地 观测 NPP 与 模拟 NPP 的 线性 拟 合 散 点 图 
Fig.2 Linear fitting scatter plot of Dangxiong 2004 to 2005 ,Haibei Meadow 2002 to 2004 ,Haibei Wetland 


2003 to 2005 overall observations and simulated values 


总 量 3.445 Gt a ^ (C) [1 16.3% 。 全 国 总 面积 2.2.2 区 域 多 年 平均 气温 、 年 降水 量 和 年 NPP 的 
为 9.59 x 10° km^ ,而 研究 区 总 的 面积 为 2.6x10' BM aN 青藏 高 原生 态 区 NPP 对 气温 变 暖和 降 
km , 占 全 国 总 面积 的 26.9% ,因此 估算 的 青藏 高 原 ” 水量 有 一 定 的 响应 ,生态 区 内 不 同 区 域 的 WPP 受气 
植被 NPP 低 于 全 国平 均 水 平 。 图 3 为 青藏 高 原 年 。 候 的 影响 而 不 同 。 

平均 气温 \ 年 均 降 水 量 和 植被 NPP 的 时 间 变 化 ,由 图 4 为 2000—2014 年 青藏 高 原 多 年 平均 气温 、 
图 3a 可 知 ,青藏 高 原 气温 起 伏 变化 较 大 ,在 2000 年 。 降水 量 和 净 初 级 生产 力 (NPP) 的 空间 分 布 格局 。 
和 2008 年 是 气温 最 低 时 期 ,但 整体 呈 略 微 增 加 的 趋 。 由 图 可 以 看 出 ,青藏 高 原 的 净 初 级 生产 力 空 间 分 布 
势 。 图 3b 显示 , 近 15 a 来 青藏 高 原 降 水 量 整体 呈 “(图 4c) 表 现 出 由 东南 部 森林 区 向 西北 沙漠 地 带 递 
下 降 趋势 ,尤其 在 2006 年 .2013 年 和 2014 年 经 历 减 的 分 布 趋势 ,该 趋势 与 研究 区 水 (图 4b) 热 (图 
了 干旱 和 较 干 旱 时 期 。 在 降水 和 气温 的 共同 作用 4a) 分布 趋势 基本 一 致 ,与 前 人 人 研究 结果 也 是 一 致 
下 , 近 15 a 来 青藏 高 原 植被 NPP 整体 呈 上 升 趋势 ， “的 后 。 在 研究 区 内 达 旺 一 察 隅 区 的 森林 区 植被 
取 值 范围 在 197.57 ~234.46 gC - m^ - a^! 其 中 NPP 最 高 在 1 000 gC .ma 以 上 ,由 于 该 区 水 
最 小 值 在 2003 年 ,最 大 值 在 2010 年 (图 3c)。 但 在 ” 热 条 件 较 好 ,气温 最 高 达 22 % ,降水 量 最 大 为 1 300 
2012—2014 年 有 明显 的 下 降 趋 势 ,为 此 ,我 们 进 一 。 mm，, 且 以 森林 植被 为 主 ,植被 覆盖 度 高 。 在 青藏 高 


步 分 析 气象 要 素 对 青藏 高 原 植 被 NPP 的 影响 。 原 东 部 和 东南 部 的 河谷 地 区 由 于 海拔 低 、 水 热 条 件 
-02;; (a) 年 均 气温 440 | (b) 年 均 降水 量 x. 240 | (QO) 年 均 NPP 
-0.4 Š 230 
只 -0.6 4 
"ET g 220 
E: -1.0 & 210 
过 -1.2 & 
> 200 
-1.4 g 
-1.6 过 190 


图 3 2000—2014 年 青藏 高 原 年 均 气 温 \ 年 均 降 水 量 和 年 均 NPP 的 时 间 变 化 


Fig.3 Time variation of annual average temperature , precipitation and vegetation NPP in the Qinghai-Tibet Plateau from 2000 to 2014 
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图 4 青藏 高 原 多 年 均 气 温 、 降 水 量 和 的 NPP 空间 分 布 


Fig.4 Spatial distribution of the average annual temperature ,precipitation and net primary productivity of the Qinghai-Tibet Plateau 


好 ,植被 NPP BEAK, TE 600 eC +m? sea AE, $ 
藏 高 原 东部 和 中 部 大 部 分 地 区 广泛 分 布 着 草原 和 草 
fa) ,这 些 地 区 的 植被 NPP 为 200 ~400 gC «m - 
a ,在 青藏 高 原 西 部 藏 北 草 旬 区 和 西南 部 灌木 区 , 
由 于 受到 水 热 条 件 的 限制 ,其 NPP 值 在 100 ~ 200 
a ”之 间 。 青 藏 高 原 藏 北部 多 为 荒漠 , 植 
被 稀少 其 NPP 值 较 小 。 
2.3 ”区域 年 均 气 温 、 年 降水 量 和 年 NPP 的 变化 趋势 
植被 净 生 产 力 NPP 的 变化 对 气候 的 响应 主要 
考虑 了 降水 量 和 气温 的 影响 。 为 了 明确 近 15 a 青 
藏 高 原 年 均 气 温 、 年 降水 量 和 年 NPP 的 变化 趋势 ， 
本 文 运用 了 Manner-Kendall ( M-K ) 非 参数 检验 法 和 
Sen's 斜率 估计 法 计算 了 青藏 高 原 年 平均 气温 .年 降 
水 量 和 年 NPP 的 变化 值 。 图 5 是 年 均 气 温 \ 年 降水 
量 和 年 NPP 的 变化 趋势 。 由 图 5a 可 知 , 近 15 a 青 
藏 高 原 大 部 分 地 区 的 气温 是 呈 增 加 趋势 的 ,尤其 是 
在 中 部 的 草 甸 区 和 东南 部 的 森林 区 ,气温 上 升 最 高 。 
图 Sb 显示 降水 量 在 中 部 的 草原 区 和 东部 的 达 旺 一 
察 隅 区 以 及 波 密 一 川西 南部 地 区 减少 ,在 青藏 高 原 
东部 及 东北 部 增加 。 在 青藏 高 原 的 西部 ,降水 量 和 


eC *m^ 


气温 都 是 减少 的 。 在 气温 和 降水 的 共同 作用 下 , 青 
藏 高 原 植被 NPP 整体 呈 增 加 趋势 ,尤其 在 青藏 高 原 
东部 及 东南 部 的 森林 区 增加 趋势 比较 大 (图 5c)。 
然而 ,在 达 旺 一 察 阳 的 森林 区 ,气温 略 有 增加 ,降水 
量 急剧 减少 ,使 得 该 区 的 植被 NPP 也 呈 减 少 趋势 。 
青藏 高 原 中 部 平原 草原 区 ,植被 NPP 随 气温 的 上 升 
而 增加 。 东 部 的 草 甸 区 ,由 于 气温 和 降水 量 共 同 增 
加 的 趋势 ,植被 VPP 也 逐渐 增加 。 青 藏 高 原 北部 多 
为 荒漠 ,植被 稀少 ,加 之 气温 不 断 升 高 ,水 分 蒸发 量 
远 远 大 于 降水 量 ,导致 柴 达 木 贫 地 区 的 植被 NPP S 
2.4 区 域 植被 NPP 与 气温 、 降 水 量 的 相关 性 

为 进一步 探讨 2000 一 2014 年 青藏 高 原 植被 
NPP 与 年 均 气 温 及 年 均 降水 量 的 关系 ,利用 Pearson 
相关 系数 法 逐 象 元 的 2 个 时 间 序 列 (NPP 与 气温 、 
NPP 与 降水 量 ) ,计算 了 近 15 a 青藏 高 原 年 NPP 与 
年 均 气温 和 年 降水 量 的 线性 相关 系数 的 空间 分 布 。 
根据 相关 系数 的 显著 性 数据 表 , 当 样本 量 为 15 时 ， 
相关 系数 的 绝对 值 大 于 0. 641 为 极 相 关 , 大 于 0. 514 
为 显著 相关 ,其 余 为 一 般 相 关 。 因 此 对 相关 系数 进 
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行 了 显著 分 类 ,如 图 6。 由 图 可 知 ,通过 显著 性 检验 
的 像 元 相对 较 少 ,上 且 相 关系 数 的 空间 分 布 与 前 人 的 
研究 基本 一 致 ” 。 

青藏 高 原 地 势 呈 西高 东 低 的 特点 。 相 对 于 高 原 
边缘 区 的 起 伏 不 平 ,高 原 内 部 反而 存在 一 个 起 伏 度 
较 低 的 区 域 。 因 其 独特 的 地 形 地 貌 使 得 高 原 区 域 气 
温 分 布 不 均匀 ,高 原 西 北部 的 羌 唐 高 原 为 极 寒 区 , 茂 
北 高 原 可 可 西里 一 带 是 高 原 的 寒冷 区 。 同 海拔 高 度 
上 来 看 ,青藏 高 原 南部 和 中 部 是 高 原 的 温暖 区 ,高 原 
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东南 部 的 气温 也 比较 高 。 图 6a 显示 , 近 15 a 来 青 
藏 高 原 地 区 基于 像 元 的 NPP 与 年 均 温 以 正 相关 为 
主 ,尤其 在 青藏 高 原 中 部 的 草原 区 表现 为 极 显 著 相 
关 , 在 东部 的 草包 区 表现 为 显著 相关 。 在 青藏 高 原 
地 区 , 正 相 关 的 地 区 占 整 个 青藏 高 原 面 积 的 
82.24% ,其 中 58. 3% 为 一 般 正 相关 ,23. 94% 为 显 
著 或 极 显著 正 相 关 。 

由 于 海拔 的 升 高 和 地 形 的 变化 ,高 原 区 域 分 布 
上 降水 量 差异 很 大 , 从 高 原 东 南部 的 森林 区 到 中 部 
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Fig.5  Qinghai-Tibet Plateau in the past 15 years annual average temperature , annual precipitation , spatial variation of NPP 
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图 6 ZIRI NPP 与 年 均 气温 、 年 降水 量 的 相关 性 分 布 


Fig.6 Correlation distribution of pixel-by-pixel NPP with annual mean temperature and annual precipitation 
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的 草原 区 ,高 寒 草 甸 区 直至 西北 的 戈壁 沙漠 地 带 , 降 
水 量 随地 形变 化 依次 减少 ,最 多 的 地 区 达 1 000 mm 
以 上 ,而 最 少 的 地 区 还 不 到 50 mm。 近 15 a 来 高 原 
降水 量 呈 现 多 元 变化 特征 ,表现 为 东部 地 区 降水 增 
加 最 明显 ,而 中 部 及 东南 部 降水 量 呈 减少 趋势 。 基 
于 象 元 的 NPP 与 降水 量 呈 正 相 关 的 地 区 主要 分 布 
在 青藏 高 原 西部 草原 区 .东部 和 南部 的 森林 区 及 西 
北部 分 草 甸 区 ,其 中 显著 或 极 显 车 相关 区 域 零星 分 
布 在 西部 .东部 及 柴 达 木 盆地 南部 。 负 相关 集中 分 
布 在 中 部 的 草原 区 . 柴 达 木 盆 地 区 和 东南 部 森林 区 , 
且 中 部 的 草原 区 和 东南 部 森林 区 表现 为 极 显 著 的 负 
相关 (图 6bb)。 青 藏 高 原 地 区 ,NPP 与 降水 量 相关 
性 大 于 零 的 地 区 占 30% ,0. 8% 为 显著 或 极 显 著 正 
相关 , 负 相 关 的 面积 占 49. 31% ,其 中 15. 95% 为 显 
著 或 极 显 著 负 相关 。 在 中 部 的 草原 区 气温 对 植被 
NPP 的 响应 大 于 降水 量 , 而 达 旺 一 察 隅 森林 区 量 对 
植被 NPP 的 影响 大 于 气温 。 总 体 而 言 ,气温 对 于 植 
被 NPP 的 影响 在 青藏 高 原 大 部 分 地 区 为 正 效 应 。 


3 讨论 


近年 来 ,植被 VPP 对 气候 变化 的 响应 是 研究 热 
点 之 一 ,不 同 的 学 者 表明 了 不 同 观 点 。 陶 波 等 ”' 通 
过 对 中 国 陆地 净 初 级 生产 力 时 空 模拟 分 析 表 明 全 国 
植被 NPP 受 厄 尔 尼 诺 现象 影响 较 大 。 姚 玉 壁 等 
研究 结果 表明 长 江 源 地 区 近 50 a di^ Ve FI" SERE 
到 “ 暖 湿 型 ”, 植 被 NPP 也 表现 为 明显 的 上 升 趋势 。 
本 文 研究 结果 表明 ,青藏 高 原 地 区 气温 逐年 增加 、 降 
水 量 逐 年 减少 的 趋势 ,同时 青藏 高 原 的 植被 年 NPP 
呈 逐 年 增加 的 趋势 ,这 与 众多 研究 结果 一 致 ” 。 

经 分 析 表 明 ,在 年 均 温和 降水 量 两 个 因素 中 , 降 
水 量 主要 影响 了 东南 部 森林 区 而 气温 是 影响 植被 
NPP 的 主要 因子 ,这 与 前 人 研究 结果 一 致 ”; E 
究 时 段 内 ,青藏 高 原单 位 面积 平均 NPP 为 217. 56 
ia ,高 于 德 吉 央 宗 等 ”的 研究 结果 。 这 
可 能 由 于 研究 时 段 不 同 ,研究 区 面积 不 同 , 且 植 被 
NPP 受 多 种 因素 综合 影响 , 除 极端 的 气候 条 件 外 ， 
还 会 受到 人 类 放牧 、 砍 伐 的 影响 ,使 得 土地 退化 , 植 
被 减少 ” 。 


eC +m 


4 结论 


本 文 基于 CEVSA 模型 估算 了 青藏 高 原 植被 净 


初级 生产 力 (NPP) 并 探讨 了 其 时 空 变化 特征 、 多 年 
变化 趋势 及 与 两 个 主要 气象 因子 的 关系 ,通过 分 析 
主要 得 出 以 下 结论 : 

(1) CEVSA 模型 通过 模拟 当 雄 、 海 北 草 乌 \ 海 
北 湿地 的 植被 NPP 进行 验证 , 当 雄 模拟 值 与 观测 值 
的 相关 性 为 0.71 ,海北 草 旬 的 相关 性 为 0. 70 ,海北 
湿地 的 相关 性 达 0. 81 ,总 体 来 看 , 3 个 站 点 模拟 效 
果 较 好 , 均 在 合理 范围 内 。 验 证 表明 ,CEVSA 模型 
在 青藏 高 原 模拟 植被 NPP 中 具有 一 定 的 适用 性 。 

(2) 时 间 分 布 上 ,气温 和 降水 量 呈 现 出 不 同 幅 
度 的 变化 ,2000 一 2014 年 间 气 温 整体 略微 增加 ,而 
降水 量 呈 下 降 趋势 。 青 藏 高 原 年 均 NPP 整体 呈 上 
升 趋势 , 取 值 范围 在 197. 57 ~ 234.46 gC - m^? - 
a ,NPP 随 气 温和 降水 量 的 变化 呈现 出 区 域 异 质 
性 ,由 东南 森林 区 向 西北 沙漠 地 带 递减 的 趋势 ,与 水 
热 分 布 趋势 基本 一 致 。 

(3) 近 15 a 来 ,青藏 高 原 西部 和 北部 气温 减 
少 ;与 气温 变化 不 同 ,青藏 高 原 中 部 、 西 部 以 及 东南 
部 地 区 的 降水 量 是 减少 的 ,在 东部 和 东北 部 明显 增 
加 。 中 部 地 区 植被 NPP 受气 温 的 影响 较 大 ,而 东南 
部 森林 区 对 降水 量 的 响应 较 大 。 

(4) 青藏 高 原 地 区 气温 ,降水 量 与 NPP 的 关系 
不 同 。 青 藏 高 原 年 VPP 与 多 年 平均 气温 之 间 呈 正 
相关 的 面积 占 整个 青藏 高 原 的 82. 24% ,年 平均 气 
温 升 高 对 植被 NPP 有 显著 的 正 向 影响 ;相反 ,NPP 
与 降水 量 呈 正 相关 的 地 区 占 30% ,年 降水 量 的 减少 
对 植被 NPP 产生 微弱 负面 影响 。 因 此 ,气温 是 决定 
NPP 增加 的 主导 因素 。 
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Response of net primary productivity of Tibetan Plateau vegetation 
to climate change based on CEVSA model 


XU Jie', CHEN Hui-ling, SHANG Sha-sha', YANG Huan’, ZHU Gao-feng , LIU Xiao-wen' 
(1 College of Earth and Environmental Science , Lanzhou University , Lanzhou 730000 , Gansu , China ; 
2 Jiangxi Provincial Water Conservancy Planning Design And Research Institute , Nanchang 330029 , Jingxi , China) 


Abstract: As the “Third Pole" of the world ,the Qinghai-Tibetan Plateau , China is located at 74 — 104 ?E and 25 - 
40 ?N. It is known for its high altitude , complex terrain ,and harsh climate. The Qinghai-Tibetan Plateau is one of the 
most sensitive regions to global climate change. Increases in temperature and changes in precipitation can affect the 
productivity of plateau’s ecosystem. Calculating the net primary productivity of the Qinghai-Tibetan Plateau ecosystem 
is very important for accurately estimating the global carbon cycle. As NEP cannot be measured directly on a regional 
or global scale ,model-based estimation is the only way to proceed. In this study ,we analyzed the temporospatial pat- 
terns and trends of NPP in the Tibetan Plateau between 2000 and 2014 using the process-based ecosystem model 
CEVSA (Carbon Exchange between Vegetation ,Soil and Atmosphere) , M-K trend test method ,the Sen’s slope esti- 
mation method ,and the Pearson correlation coefficient method. The model was based on a 0. 1? x0. 1? resolution map 
of vegetation types ,soil texture data , and daily meteorological data. The results indicated as follows: (1) the net prima- 
ry productivity for the Qinghai-Tibetan Plateau decreased from southeast to northwest and was consistent with the trend 
of water-heat distribution. The results were similar to those obtained by ZHOU Caiping et al, who applied a combina- 
tion of terrestrial ecosystem model and MODIS data to estimate the net primary productivity of the Qinghai-Tibetan 
Plateau. In the spatial distribution ,the NPP of the forests in the east and southeast was between 600 and 1200 gC - 


2 


m ^-*a'!,the NPP of the central grassland and meadows was from 200 to 400 gC - m^ 


2 


- a !,and in the western 
and northern deserts ,the NPP was limited by the moisture and temperature. (2) the annual average temperature in- 
crease had a significant positive effect , while the precipitation decrease had a significant negative effect on the Qing- 
hai-Tibetan Plateau’s NPP. The annual NPP was positively correlated with annual mean temperature over 82.24% of 
the region , while negatively correlated with annual precipitation over 49.31% of the region. Therefore , temperature is 
considered to be the dominant factor determining spatial variations in NPP. The predecessors also obtained similar re- 
sult. For example , LIU Gang et al. analyzed the spatiotemporal variation of net primary productivity and climate con- 
trols in China from 2001 to 2014. Based on their results ,the correlation analysis between NPP and meteorological ele- 
ments indicated that NPP was positively correlated with temperature in the Changbai Mountain area , Qinghai- Tibetan 
Plateau ,and southern areas. (3) from 2000 to 2014, the trend of increasing NPP was consistent with the changes in 
temperature. The precipitation showed a slight decrease change. A period of warming accompanied by a decrease in 
precipitation contributed to the trend of a gradual increase of NPP in Qinghai-Tibetan Plateau. Therefore , improving 
our ability to accurately describe the response of NPP to climate changes will provide a better understanding of terres- 
trial ecosystem responses to global changes. However , there are still some issues that need to be addressed ,such as the 
uncertainty of NPP prediction. These uncertainties mainly include the driving variables and parameters in the model. 
Overall ,in order to predict the impact of climate change on ecosystems at the regional level , modeling uncertainty can 
be reduced by increasing the spatial resolution of the driving variables. In addition, optimizing the model parameters 
can also reduce the uncertainty in the model simulation. 


Key words: CEVSA model; net primary productivity; climate change; Tibetan Plateau 


